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Abstract 

The formation of bulk metallic glass requires the constituent elements to have a negative heat of mixing 
but has no restrictions on its magnitude. An understanding of this issue is lacking due to the absence of 
a valid method for describing chemical ordering of metallic glasses. For example, the radial distribution 
function is ineffective in identifying the elemental preferences of packed atoms. Here, we show that using 
molecular-dynamics simulation, the chemical medium-range ordering of liquid alloys can be evaluated from 
persistent homology. This inherently arising chemical medium-range order in metallic glasses is exclusively 
regulated by the activation and inhibition of the constituent components, making the topology of metallic 
glasses a Turing pattern. The connecting schemes of atoms of the same species form three distinct regions, 
reflecting different correlations at the short and medium length scales, while the difference in the schemes 
corresponds to chemical ordering. By changing the elemental types, it is demonstrated that the chemical 
medium-range order strongly depends on the relative depth of the interatomic-potential wells. The study 
separates metallic glasses from crystals under the condition of negative heat of mixing by emphasizing their 


fundamental difference in interatomic potentials. 


INTRODUCTION 


According to the empirical rule of Inoue, one of the prerequisites for casting bulk metallic 
glasses (MGs) is the negative heat of mixing (AHmix) between any two elements of the con- 
stituents [1]. A negative A Hmix suggests a bonding choice preferential to unlike species against 
like species, so factors like phase separation that are detrimental to glass formation during casting 
can be avoided. One might infer from this that giant MGs can be produced by lowering AHmix 
continuously or by choosing elements with A Hmix of higher absolute values. Unfortunately, the 
truth is different: a study of 7688 alloys revealed that the critical cooling rate (R,) of binary MGs 
is only weakly correlated with the absolute value of negative AHmix [2, 3]. It is unclear why a 


negative A Hmix is significant but its magnitude has no bearing on MG formation. 


Since mixing involves chemical affinity of the elements, the task requires an understanding of 
the chemical ordering inside the MGs [4, 5], which implies bonding preference and how far the lo- 
cal chemistry departs from a random melt. Our knowledge, however, is rather limited and largely 
focused on the first radial-distribution-function (RDF) peak. Although many ordering characteris- 
tics describing the structural and chemical order of the first nearest neighbors are proposed [6-14], 
it is unclear what ordering characteristics contribute to the RDF peaks after the first one. Addition- 
ally, although the chemical affinity of the solute and solvent atoms of alloys are resembled using a 
variety of radii by Miracle [7, 15-17] following up with Bernal’s hard sphere model [18-22], these 
approaches are still insufficient to resolve the mystery of negative AH,,;, in the formation of MGs. 
A model emphasizing the chemical ordering is required, because configurations would change if 
the affinity of the alloying elements changed even though the radii are set. To understand the 
medium-range ordering (MRO) of amorphous materials, numerous simulation studies have been 
carried out [23-25], but they all have some drawbacks. The bond-angle or dihedral-angle distri- 
bution only highlights configurations of the third-nearest neighbors [26, 27]. Ring statistics are 
only useful for polymeric glasses and crystals [28—30]. The MRO topology was defined in earlier 
investigations of persistent homology analysis, but the chemical influence was missed [31, 32]. On 
the other hand, an analogy from DNA helps us to understand the significance of finding CMRO in 
liquids. While the statistical findings from X-ray diffraction assisted people in understanding the 
double helix structure of DNA reflecting MRO, the sequence of complementary base-pairing of 
DNA representing chemical order, lays the foundation for the variation of life. Thus, developing a 


new model that visualize the chemical medium-range order (CMRO) is what waits to be resolved. 


2 


The current work provides a simple mathematical model that, while ignoring the specifics, 
draws out the logical connections in complex glass-forming liquids. The concept is derived from 
Turing’s reaction-diffusion theory [33, 34], which explains how ordered patterns spontaneously de- 
velop from random complex systems with only an activator and an inhibitor and has been applied 
in many other fields, such as biology [35-37], chemistry [38—41], and 2D materials [42, 43]. In de- 
tails, the elemental distribution of Cus oZrsq metallic glass is firstly visualized through continuous 
patterning, which is then confirmed as a Turing pattern. Secondly, persistent homology analysis 
is applied to study how temperature and interatomic potential influence the topology. Thirdly, a 
new criterion in addition to the negative AH mix is proposed, in order to separate the conditions for 
making metallic glasses from making crystals. Persistent homology analysis demonstrates that the 
relative magnitudes of the well depth of the interatomic potentials determines the medium-range 


structure of liquid alloys, which separates the products into categories including metallic glasses. 


RESULTS 


Modeling and its verification 


FIG. 1: The stucture of a 2D Cuso Zrso liquid. (a) Red (for Cu) and blue (for Zr) balls 
highlighting the atomic positions. (b) Red (for Cu-Cu) and blue (for Zr-Zr) lines highlighting 
atomic bonds. (c) Continuous patterning highlighting the chemical ordering. The dimensions of 


the representation are 5 x 5 nm/?. 


A reformative microscopic depiction of metallic melts may help us grasp the CMRO. CusoZrs0, 


a model MG-forming liquid, is used as an example. The model is a 45.6 x 45.6 nm? 2D square 
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consisting of 29600 atoms. The sample was equilibrated at 2000 K after 1 ns relaxation. A portion 
of the square, i.e., 5 nm x 5 nm, is cut out and shown in Fig. 1. Discrete dots emphasizing atomic 
positions of Cu (red dots) and Zr (blue dots) is presented in Fig. la. However, this presentation 
has the flaw of ignoring the atomic interactions that should take up the white space in Fig. la. As 
an alternative approach, the bondings of Cu-Cu (red lines) and Zr-Zr (blue lines) atomic pairs are 
demonstrated as the colored lines in Fig. 1b, from which different morphologies of Cu and Zr are 
revealed. Here, the Cu areas are separated, whereas the Zr areas are interconnected. Because Cu 
and Zr atoms have an equivalent atomic proportion of 50%, such morphological differences are 
inherent in chemistry. Further, by filling the white space of Fig. la using cubic interpolation (see 
more details in Fig. S1) , the Cu and Zr areas are highlighted in the manner of Fig. 1c emphasizing 
the morphological difference between Cu and Zr. In this depiction, the minutiae of atomic position, 
atomic size and bonding are ignored, and the gist of chemical distribution is highlighted. These 
morphological traits are referred to as the chemical ordering of the Cus pZrso liquid. 

Figure 1c resembles a Turing pattern. To determine a Turing pattern, however, four require- 
ments must be met: (i) the system consists of an activator that provides short-range positive feed- 
back and an inhibitor that provides long-range negative feedback [34, 42, 44, 45]; (ii) the static 
image is obtained from dynamic evolution [36]; (iii) different initial states have no effects on the 
topological properties of the final structure [34, 35]; and (iv) the correlation length is fixed and 
independent on sample size [38, 39]. We now look into the validity of four criteria in the CuspZrso 
liquid. 

The activator and inhibitor are determined by the interatomic potentials. Depending on the 
valley depth of the potential functions (| £z,-z,| for Zr-Zr bond, |Ecu-zr| for Cu-Zr bond and | Ecu-cu| 
for Cu-Cu bond), the chemical affinity between atoms can be evaluated. For Cu-Zr, potential depth 
decreases in the order of Zr-Zr, Cu-Zr, and Cu-Cu. The formation of Zr-Zr and Cu-Zr bonds are 
active, promoting self-activation of Zr-Zr bonds and explaining the inter-connected morphology of 
Zr (Fig. 1c). The passive formation of Cu-Cu bonds, on the other hand, results in the long-range 
stimulation of Cu and the formation of Cu-rich regions. However, because the Cu-Zr bonds are 
energetically favored, meaning the creation of two Cu-Zr bonds at the cost of breaking Cu-Cu 
and Zr-Zr bonds are exothermic (A Hix = —23 kJ mol! [1]), the Zr region is prevented from 
growing bigger, while the Cu regions are separated.In other words, Cu-Zr bonds inhibit Zr areas 
from growing over long distances, restricting Zr-Zr activation to short distances. Here, AH mix < 0 


functions in the prevention of Zr-Zr region enrichment. Fig. 2 summarizes the activation and 
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FIG. 2: The formation mechanism of Turing pattern. Here, the activator and inhibitor are Zr 


(blue balls) and Cu (red balls), respectively. 


inhibition schemes of Cu and Zr. 

A Turing pattern should be independent on the initial distribution of the activators and in- 
hibitors. Otherwise, diffusion of the activator can create a pseudo-Turing pattern, driven by the 
chemical gradient of a similar pattern. The essence of a Turing pattern is the reaction between the 
activator and inhibitor that can erase the memory of the initial morphology. Starting from differ- 
ent initial distributions of CusoZrso, we found that the stabilized morphology is always a Turing 
pattern (Fig. 3). Though the final patterns are not exactly the same, the morphology of the blue 
and red areas are similar. Movies S1 and S2 also prove that the static Turing pattern of CuspZrso is 
obtained by dynamic evolution. Further, the size of the model system does not matter with the pat- 
tern. In another sample with dimension of 22.7 x 22.7 nm? consisting of 7360 atoms, the resulting 
pattern of Movie S2 and Fig. S3 are similar to those shown in Fig. 3. Thus, the morphology of 


CusoZrso liquid meets all the requirements of a Turing pattern. 


Pattern analysis through persistent homology: temperature effects 


The above-mentioned results of 2D Cus Zrsq liquid provide fundamentals of a Turing pattern. 
We now move on to a 3D CusgZrsy liquid for topological analysis. A model system in a cube 
form of 17 x 17 x 17 nm’, consisting of 256000 atoms, is simulated. A cube with a side length 


of 2 nm is cut out, highlighting the morphologies of Cu (red) and Zr (blue) as shown in Figs. 4a 
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a Initial state I b Initial state IT C Initial state TT 


d Final state I e Final state II f Final state II 


FIG. 3: Morphology evolution of Cus 9Zrso starting from three initial conditions including 
(a) state I: Cu (red) and Zr (blue) are completely separated; (b) state II: Cu and Zr are arranged 
alternately in stripes; (c) state III: Cu and Zr are randomly arranged. The final states are presented 
in (d), (e) and (f), corresponding to (a), (b) and (c) after annealing the samples at 2000 K for (d) 
20 ns, (e) 1 ns and (f) 1 ns. Yellow double-ended arrows visualize the characteristic lengths of the 


chemical medium-range order. The dimensions of the representation are 5 x 5 nm’. 


and 4b. Here, Cu-Cu bonds are drawn when the interatomic spacing between Cu and Cu is within 
the cut-off distance (deut), and deut is set as the locations of the first valleys in RDF (as shown in 
Fig. S2 and as collected in Table. S1). The Zr-Zr bonds are drawn in the same way. For better 
presentation, the bond width, i.e., the cylinder diameter, is set at 2 A. A slice of 5 -nm-long, 5-nm- 
wide and 5-A-thick is cut out of the 3D cube, and the resulting patterns (Figs. 4c and 4d) resemble 
characteristics of 2D Turing patterns (Fig. 1c). A quantitative investigation of CMRO on a Turing 
pattern, particularly in 3D [46], however, necessitates persistent homology. Persistent homology 
is a topological method for extracting detailed properties in 3D structures [47]. It has been applied 
to extract not only structural short-range order but also structural medium-range order of oxide 
and metallic glasses [31, 32]. Unlike a former analysis [31] that neglected the element types and 
focused on the topology atomic structure for persistent homology analysis, the models applied 
here consider Cu and Zr atoms separately (Figs. 4e and 4f). 

A persistent diagram (PD) is a birth-death time diagram obtained from persistent homology 


analysis. The birth time (tp) corresponds to the interatomic distance that forms a closed ring of 
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FIG. 4: Structure analysis of CusoZrso in 3D at (-1) 3000 K, (-2) 2000 K and (-3) 300 K, 
respectively. 3D structures of (a) Cu-Cu and (b) Zr-Zr bonds in CusoZrso (2 x 2 x 2 nm). Slice 
morphologies of (c) Cu-Cu and (d) Zr-Zr bonds (5 x 5 x 0.5 nm?). Persistent diagrams of (e) Cu 

and (f) Zr. ta, tp, Lo and deut denote the death time, birth time, correlation length and cut-off 


distance, respectivelly. 


bonds, while the death time (tą) is the interatomic distance that completely divides the ring by 
triangles (see Methods for more details). Three characteristics regions are found on the PDs of 
Cu. The S-region is located at (2.9 A, 3.2 A) representing (tp, ta), both are smaller than deut of 
3.5 A, indicating that the S-region represents the rings formed by the first-nearest neighboring Cu 
atoms. Here, deut is the abscissa of the trough of the first neighbor atom (see Fig. S2 and Tab. S1 
for more details.). The M-region shown at (4.4 A, 4.8 A) represents the rings that are formed by the 
second-nearest neighboring Cu atoms. Examples are provided by “circle” Cym in Fig. 4c-1. The S- 
and M- regions have their t, close to tg. The I-region centered at (3.2 A, 5.1 A), reflects the rings 
that are formed by the first-nearest but terminated by the second-nearest neighboring Cu atoms. 


Examples are given by “circle” Cy in Fig. 4c-1, whose tg is 59% greater than its tp. Because it 
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includes a ring with tg > deut, both I and M regions contain medium-range information. Fig. 4c-1 
provides examples of S-, M- and I-regions. The S-region reflects rings formed by the first-nearest 
neighboring atoms at shorter length scale; the M-region is contributed by the rings established 
by connecting two separated Cu regions; the I-region represents rings formed by the first-nearest 
neighboring atoms at longer length scale. For the PD of Zr (Fig. 4f-1), two distinct zones are 
found. Within the deut of 4.2 A, the S-region centers at (3.4 A, 3.7 A), reflecting short-range order. 
The I-region centers at (3.3 A, 4.8 A), reflecting medium-range order. The absence of an M-region 
distinguishes the PD of Zr from the PD of Cu. This indicates there are no rings produced by two 
separated Zr regions (Fig. 4d-1). This makes sense because all the Zr atoms are linked together. 
As a result, the existence of M-region established the distinction between Cu and Zr, reflecting 
chemical ordering at the medium range. Since no topological regions can be found beyond 5.1 
A (the maximum value of coordinates of all topological regions), the CMRO as revealed by the 
persistent homology limits to the second-nearest neighboring atomic distance. The mean tq for I 
and M regions was used as the Lc of CMRO. Here, the Lo that defines CMRO is similar to both 
Cu and Zr atoms, of 5.0 A on average. The width of the stripes or the spacing between the stripes 


of Turing pattern of MG (Fig. Ic) visualizes the Lc. 


The temperature effects are further investigated. The Turing patterns exhibit similar features 
in both 3D (Figs. 4a and 4b) and 2D (Figs. 4c and 4d) at all three characteristic temperatures. 
However, the central colors and positions of the characteristic regions of their PDs differ. The 
central color of the regions warms up when the temperature is reduced from 3000 K to 2000 K, 
indicating that the liquid structure is more heterogeneous at the lower temperature. Additionally, 
the centers of all the regions shifts to the bottom left, indicating a compact local structure. At 300 
K, the influence is increasingly pronounced as the regions highlighted by the dashed curves are 
nearly all filled by the central color, and the central positions are further reduced. Additionally, 
the boundaries separating S-, I- and M- regions are more defined in 300 K (Figs. 4e-3 and 4f-3) 
compared to 2000 K (Figs. 4e-2 and 4f-2). The Lo of CMRO, on the other hand, is slightly reduced 
to 4.7 A at 300 K, in comparison with 4.8 A at 2000 K and 5.0 A at 3000 K. The structural ordering 
brought on by vitrification is reflected in a modest shortening in Lc and a narrower distribution of 


the characteristic local regions. 
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FIG. 5: Structure analysis of different binary metallic glasses in 3D at 300 K. Radar map 
characterizing the potential well of (a) Cu-Ag, (b) Cu-Zr, (c) Ni-Zr and (d) Ni-Al. Eu and r, are 
normalized well energy and distance, respectively. Persistent diagrams and slice structures 
(5 x 5 x 0.5 nm?) of (e) CusgAgso, (£) CusoZrso, (g) NisoZrso and Nis Also. Here, red, gray, blue, 
yellow and green represent Cu, Ag, Zr, Ni and Al respectively. The upper-left triangle is the 


persistent diagram, and the bottom-right triangle is the slice structures. 


Pattern analysis through persistent homology: Elemental effects 


The effects of elemental types are studied on four equiatomic compositions, i.e., CuspAgso, 
Cus 9Zrs0, NisoZrso and NisoAlso. Their radar maps (Figs. 5a, 5b, 5c and 5d) represent variations 
in the normalized depth of potential well (Eu = ne) and the normalized distance of potential 


well (ru = ). Here, E and r respectively are the valley energy and distance of the potential 


Tmax 


function, while Emin and Tmax respectively denote the lowest values of E reflecting the deepest 
valley and the highest r among all the alloy constituents. So, a stronger atomic binding is denoted 
by a larger Fu, and a shorter atomic bond is denoted by a smaller rą. For CuspAgso (Fig. 5a) and 
CusoZrso (Fig. 5b), the rą-triangles are comparable, but the /,-triangles are different. The Cu-Cu 
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bond is the strongest in Cusy,Ag59 and the weakest in CusoZrsoọ. Due to this, the S-region in the PD 
of Cu in CuspAgso (Fig. 5e) and the M-region in the PD of Cu in CusoZrso (Fig. 5f) are both the 
strongest. However, the PD of Ag in CusgAgso (Fig. 5e’) resembles the PD of Zr in CusoZrso (Fig. 
5f’), both containing S- and I-regions. Despite having various E,-triangles, the similar topological 
properties are determined by the large potential well distance of Ag and Zr (7,=1). Given that the 
S-regions of Cu and Ag are both the warmest in color among their S-, I- and M-regions, it suggests 
that both the activator and inhibitor of CuspAgso are closely packed in topology (Figs. 5e and Se’), 
forming a microscopic structure of phase separation, and creates difference from the topology of 
CusoZrso (Figs. 5f and 5f’). This is expected by the positive AH,,;, of Cu-Ag[1] and supported by 
experimental evidence [48, 49]. In the case of NisgZrs, it likewise possesses a similar r,,-triangle 
to CusoZrs0, but the Ni-Zr bond is stronger than Cu-Zr bond in terms of Eu (Figs. 5b and 5c). As a 
result, the I-region is stronger in the PDs of both Ni and Zr (Figs. 5g and 5g’) than it is in the PDs 
of both Cu and Zr (Figs. 5f and 5f’). The contradiction comes from Ni-Al which has a negative 
AHnmnix [1], and is a bad MG former as reported[50]. The Ni-Al bond is the strongest for NisoAlso 
(Fig. 5d), therefore it tends to break Ni-Ni and Al-Al bonds. This rises the intensity of region I for 
both Ni and Al (Figs. 5h and 5h’), resulting in a binary crystal structure easily. Thus, the clues of 


glass formation resides in the intensity order of these PD regions. 


DISCUSSION 


We now provide a schematic explanation (Fig. 6) on how chemical affinity influences the 
solid phase by using our data of equi-atomic binary alloys (As59Bs50). These alloys have similar 
ry-triangle but different combinations of interatomic potentials. Here, the chemical affinity is 
represented by the bond energy or the depth of the interatomic potentials (€ < 0). The absolute 
value of E of A-A bond, | £44], is assumed greater than that of B-B, i.e. |E4a| > |Eppl. 

If AHmix > 0 occurs when |Eaa| > |Egg| > |Fap|, full immiscibility results in the sample. 
If AHmix > 0 is obtained at |EĘaa| > |Eap| > |Epp|, the system is partially miscible, and phase 
separation occurs. Here, B-B bonding is the least favored so that B-rich regions are hindered. 
On the other hand, A-A maintains its chemical affinity to bond with each other, which allows for 
A-A bond enrichment. Because A-B bond is stronger than the B-B bond, the A-rich regions are 
separated by the A-B bonds and surrounded by the passively formed B-rich regions, leading to 


phase separation. In the case of AHmix < 0, if it is acquired by |Eaa| > |Eap| > |Epp|, MGs 


10 


Energy Product Intensity Order 
Element: A B 
Immiscible crystals SLM, Sp>l5>Mp 


FIG. 6: Schematic diagram explaining how interatomic-potential energy depth (EAA, EBB 

and Fas) and the intensity order of PD regions influence the products. Because the heat of 

mixing is AHmix = EAB — aatia, 51], the products are categorized into four sections,as 
separated by the dashed lines, according to the relative position of Eag. Sa (Spg), Ia (Ig) and 


Ma(Msg) represent the peak intensity of PD regions S, I and M of element A (or B), respectively. 


are formed. The A-A bonding are still preferred as A-A bond has the strongest chemical affinity, 
but the growth of A-rich regions are strongly inhibited due to the growing influence of A-B bonds, 
making crystallization from A-rich regions difficult. If, on the other hand, AH mix < 0 is obtained 
at |E'ap| > |Eaa| > |Epp|, A-B bonding is energetically preferred so that A-B bond rich regions 
are grown. Crystalline binary alloys are likely produced at this condition. 

The PD analysis provides a structural basis for these observations. When A and B atoms are 
fully immiscible, they are separatedly packed, and S4 and Sz are the strongest on their PDs. For 
phase separation, Sa and Spg are still the strongest because of the active enrichment of A atoms 
and the passive enrichment of B atoms. But, because the B-rich regions are passively formed, 
lacking connections to each other, a strong Mpg resembling linkage of isolated regions in topology 
is obtained, giving rise to the relationship Sg > Mg > Ip. In the case of MG formation, A atoms 
are activators whose intensity order remains unchanged with S, being the strongest. However, 
since B atoms are bonded more strongly with A atoms, the B-rich regions are capped by A atoms, 
leading to a spatial distribution of discrete islands. This then makes the Mg the strongest in its 


PD. Lastly for crystalline alloys, I, and Ip are the strongest on their respective PDs, reflecting 


11 


the influence of strong A-B bonding. In this case, although AHmix < 0 its magnitude does 
not promote MG formation but instead drives crystallization. Therefore, it is crucial to maintain 
|Exa| > |Eap| > |Esg| while keeping a negative AHmix to create MGs. This is the criteria 
resolved from our topological analysis. 

Finally, we highlight the broad implication this research may have. The basis of the model is 
solely dependent on the binding-energy order, regardless of the kind of objects or the size of the 
system. The fundamental building blocks can be applied to any objects like clusters or celestial 
entities in addition to atoms. In theory, any systems with competitive relations might have topo- 
logical patterns in common. The different ways that the constituents engage with one another are 


what bring order out of chaos in spatial arrangement. 


CONCLUSION 


This paper develops a new model to visualize the chemical medium range order of metallic 
glasses and answers why the magnitude of negative heat of mixing is irrelevant for glass forma- 
tion. By shifting the microscopic presentation from discrete to continuous patterns, the present 
work discloses a Turing pattern with traits of dynamic development, starting-structure irrelevance, 
and sample-size independence. The patterns are quantitatively analyzed by persistent homology, 
where the topological difference of the constituents are revealed. Temperature and elemental type 
influence the topological properties shown in the persistent diagrams. Apart from a negative heat 
of mixing, the bonding energy of unlike species has to be set in the middle of the bonding en- 
ergy of like species so that the intensity order of the persistent-homology regions is in favor of 


metallic-glass formation. 


METHODS 
MD simulations 


Molecular dynamics (MD) simulation was conducted on LAMMPS [52] software. In this pa- 
per, the potential functions of Cu-Zr [53], Cu-Ag [54], Ni-Zr [55] and Ni-Al [56] were used. An 
isobaric and isothermal (NPT) ensemble was applied, and the system pressure was set to 0 with 
periodic boundary conditions. 2D CusoZrso, containing 29600 atoms in a 45.6 nm x 45.6 nm 


square, and all 3D samples, containing 256000 atoms in cubes, were modeled. The 2D and 3D 
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liquids were simulated by annealing the specimens at 2000 K and 3000 K for 1 ns to ensure that 
their liquid configurations were uniform. The 3D glass were simulated by quenching the liquids 


from 3000 K to 300 K at a cooling rate of 2.7 x 10! K s™t. 


Persistent Homology 


Persistent diagrams were constructed to analyze the topologies of different kind of atoms. The 
principles [31, 32] were to determine the birth time (¢;,) and death time (tą) of virtual rings. Due 
to the software’s restriction on the number of atoms, 5000 atoms, carved from the center of the 
3D specimen, were employed. All the persistent diagrams in this paper were obtained from the 
analysis of about 5000 atoms. In this work, Ripser toolkit [57] based on the Vietoris-Rips [58] 
was employed. In this toolkit, the virtual rings were born by connecting the virtual circles orbiting 
the atoms and were died when entirely divided by triangular loops. An example of t» and ta 


determination was given in Fig. S4. 
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S1 


SUPPORTING MOVIES 


Movie S1 shows the pattern evolution of a 2D 29600-atom Cus Zrs9 system at 2000 K with an 
initial state Cu and Zr are completely separated. 
Movie S2 shows the pattern evolution of a 2D 7360-atom CusoZrs9 system at 2000 K, with an 


initial state when Cu and Zr are sliced into multilayers. 


SUPPORTING INFORMATION 


A grid meshes the discrete representation of atomic positions. When the center of a Cu atom 
falls within a grid element, the grid element is indicated as 0, and when the center of a Zr atom 
falls within a grid element, it is indicated as 1. The remaining grid elements are filled by “cubic” 
interpolation of gridata code in Matlab 2018b. After that, the numerical matrix is binarized by 
setting the entries to 1 or 0 depending on whether the value is equal to and larger than or smaller 


than 0.5. As a result, a continuous pattern is created, with red denoting elements at values of 0 and 


blue denoting those at values of 1. 
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FIG. S1: Schematic diagrams of the continuous process of 2D Cu50oZr5o structure. 
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Figure. S2 shows the g(r) of CusoZrso at different temperatures in 2D and 3D. 
The positions of the 1* peak and its following valley are summarized in Table S1, where the 


positions of the valley defines the cut-off distance (deut). 


Time (t) is visualized by the diameters of the circle centering each atom. When two circles 
are touched, a virtual bond is establish with the bond-length equal to the diameters. In this work, 


Ripser toolkit{S1] based on the Vietoris-Rips[S2] was used. When the virtual bond is connected 
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FIG. S2: Radial distribution function (RDF) g(r) of 2D CusoZrso at (a) 2000 K, (b) 300 K and 
3D CusoZťrso at (c) 2000 K, (d) 300 K. 


Table S1: The locations of peaks and valleys in RDF under different circumstances. 


State | Cu-Cu (Å) | Zr-Zr (A) | Cu-Zr (A) 


2D peak|valley | peak | valley | peak | valley 


300K | 2.5 | 3.4 | 3.0 | 4.0 | 2.6 | 3.8 


2000 K| 2.4 | 3.5 | 3.0) 4.1 12.6] 3.8 


3D peak|valley | peak | valley | peak | valley 


300 K | 2.6 | 3.3 | 3.2 | 4.0 | 2.8 | 3.8 


2000 K| 2.5 | 3.5 |32 | 4.2 |27| 3.8 


3000 K| 2.5 | 3.5 | 3.2 | 4.2 |2.7| 3.8 


into a polygon (quadrilateral and above), a ring is born, this time is recorded as birth time (tp). 
When the polygon is completely split into a combination full of triangles, the ring is dead, and 
this time is denoted as death time (tq). For example, as shown in Fig. S4a, the C: ring consisting 


of four neighboring bonds is created at t = dı, yet ended at t = dz because its ring structure is 


S3 


FIG. S3: Turing structure in a 7360-atom Cus 9Zrs59 system after annealing it at 2000 K for 1 ns, 
which is consistent with that of 29600-atom system after annealing it at 2000 K (Fig. 3), 


indicating that the size of the system does not affect the morphology. 


completely split into triangular loops. 

Figure S4b is called persistent diagram (PD), which takes t, and tg obtained in Fig. S4a as 
abscissa and ordinate, respectively. In Fig. S4b, a time t = deut is selected, which divides the 
PDs into 3 dashed regions: blue, yellow and red, so called region D (death), region B (birth) and 
region F (future) respectively. When t = deut, as shown in Fig. S4a, the dead ring appears in 
region D, the undead ring appears in region B, and the ring about to be born appears in region F. 
For example, when d3 < deut < d4, the connection status of the virtual bond is the same as that of 
t = d3. At this time, in Fig. S4a, Cı has been dead, Ca has been born but undead, and Cg has not 
yet been born. t» and ta of C1, C2 and Cg are drawn as blue, yellow and red points, respectively, 


in Fig. S4b, which fill into region D, B and F, respectively. 


S4 


Ss ——— 


Death time 
A 


N 


FIG. S4: Schematic diagram of persistent homology. (a) The schematic explaining the birth 
and death of the rings. Here, d is the interatomic distance. (b) Persistent diagram of (a). The blue, 


yellow and red area represent the death, birth and future regions, respectively. 
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